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Abstract 

We investigate the ergodicity of 2D large scale quasigeostrophic 
flows under random wind forcing. We show that the quasigeostrophic 
flows are ergodic under suitable conditions on the random forcing and 
on the fluid domain, and under no restrictions on viscosity, Ekman 
constant or Coriolis parameter. When these conditions are satisfied, 
then for any observable of the quasigeostrophic flows, its time average 
approximates the statistical ensemble average, as long as the time 
interval is sufficiently long. 
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1 Introduction 



The models for geophysical flows are usually very complicated. Simplified 
models have been developed to investigate the basic key features of large scale 
phenomena. These models filter out undesired high frequency oscillations in 
geophysical flows and are derived at asymptotically high rotation rate or 
small Rossby number. 

An important example of such a geophysical flow model is the quasi- 
geostrophic flow model ITS 



Aip t + J(ip, Aip) + (3ip x = vA 2r ip — rAvp + wind forcing , 

where ifj(x,y,t) is the stream function, (3 > is the meridional gradient of 
the Coriolis parameter, v > is the viscous dissipation constant, r > is 
the Ekman dissipation constant. Moreover, J(f,g) = f x g y — f y g x denotes 
the Jacobian operator. 

The quasigeostrophic equation has been derived as an approximation of 
the rotating shallow water equations by the conventional asymptotic expan- 



sion in small Rossby number fIB |. Recently, the randomly forced quasi- 
geostrophic flow model has been used to study various phenomena in geo- 
physical flows under uncertain wind forcing [|TJ|, [16|, |20|, [14], |5[ . 

Introducing (relative) vorticity u(x, y, t) = Ai/j(x, y, t), the quasigeostrophic 
equation can be written as 

uj t + J(ip, oS) + (3ip x = vAuj — ruj + wind forcing , (1) 

where (x, y) G D and D C R 2 denotes a bounded domain with sufficiently 
regular boundary. Potential vorticity is defined as u> + (3y. The boundary 
conditions are no normal flow (if) = 0) and free-slip (uj = 0) on dD as in 
Pedlosky ([|19|, p. 34) or in Dymnikov and Kazantsev ||: 



iP = lu = ondD . (2) 

An appropriate initial condition uj(0) is also imposed. We note that the 
Poincare inequality holds with these boundary conditions. 

An invariant measure for stochastic systems is like a "statistical steady 
state" and is a part of the asymptotic permanent regime of the system 
When there is only one invariant measure for the quasigeostrophic flows mod- 
eled by @, we have the so-called ergodic principle, i.e., for any observable 
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of the quasigeostrophic flows, its time average on [0, T] approaches the sta- 
tistical ensemble average, as T goes to infinity. 

We will investigate the existence and uniqueness of invariant measures 
for quasigeostrophic flows. After review mathematical setup in §2, we study 
exixtence and uniqueness of invariant measures in §3 and §4, respectively. 
Finally, we summarize our results in §5. 

2 Mathematical Setup 

In the following we use the abbreviations H = L 2 (D), H$ = Hq(D), H k = 
H k (D), < k < oo, for the standard Sobolev spaces. Let < ■, ■ > and || • || 
denote the standard scalar product and norm in L 2 , respectively. Moreover, 
the norms for Hq are denoted by || • \\H k - Due to the Poincare inequality ||, 
||A(£>|| is an equivalent norm for Hq. It is well-known that the linear operator 

A = v A : H -> H 

with domain D(A) = H 2 C\Hq is self-adjoint. Note that A generates a strongly 
continuous, and in fact, an analytic semigroup S(t) on L 2 ([0). The spec- 
trum of A consists of eigenvalues > Ai > A 2 > A 3 > . . . with corresponding 
normalized eigenf unctions e\, e2, .... The set of these eigenfunctions is com- 
plete in L 2 . For example, for the square domain D = (0, 1) x (0, 1) the 
eigenvalues are given by —vim 1 + n 2 )n 2 for positive integers m,n, and the 
associated eigenfunctions are suitable multiples of sin(m7rz) sin(n7r?/). 
We define the nonlinear operator F by 

F(lo) = -ru - Pi> x - J(V>,w), 

then (JTJ) can be rewritten as the abstract evolution equation together with 
initial condition 

du = {Auj + F{u))dt+ ^/QdW , (3) 
o>(0) given, (4) 

where W(x, y, t) is a Wiener process defined on a probability space (Q, J 7 , P). 
The covariance operator Q : H — > H for this Wiener process is a nonnegative 
and symmetric linear continuous operator to be specified below. The term 
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with Ito derivative, y/QdW , is a model for the white-in-time noise represent- 
ing the random wind forcing. This equation can be rewritten in the mild 
(integral) form 

u(t) = S(t)u(0) + I S(t- s)F(u(s))ds + Z(t) , (5) 
Jo 

where Z(t) is the stochastic convolution 

Z{t)= [ S(t- s)y/QdW(s) , t>0. (6) 
Jo 

In fact, Z(t) is an Ornstein-Uhlenbeck process and it is the solution of the 
linearized version of the above equation (|3|): 

dZ = AZdt+^QdW. (7) 

In this paper, we always assume that the covariance operator Q for the 
Wiener process W(t) is of trace class, i.e., Trace Q < +oo. Thus we only 
consider noise that is white in time but colored in space. Then the stochastic 
convolution Z(t) has a continuous version with values in H = L 2 (D); see 
Theorem 5.14 in 

We can specifically define an appropriate class of Wiener processes W(t) 
satisfying the above condition. Let for positive integer k, denote a 

family of independent real-valued Brownian motions. Furthermore, choose 
positive constant such that 

oo 2 

Er^f-<°° ( 8 ) 

for some < 7 < 1. Then we define the white noise by 



y/QW(t) := <*kPk(t)e k , t > . (9) 

k=l 

Note that the eigenvalues for the operator A behave like k in two di- 
mensions and also note that the Riemann zeta function £(s) = Ylk=i F 
is well-defined for s > 1. We see that the condition (ffl) is satisfied when 

13 

k~2 < at < A;~s. 
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We further assume that 

K (D) = iuf iuf rneas(DnB(x,y;p)) 

0<p<diam(D) (x,y)eD p 2 

where diam(D) is the diameter of D (the least upper bound of two-point 
distances in D), meas(-) denotes the Lebesgue measure, and B(x, y; p) is the 
open disk centered at (x, y) and with radius p. We also assume that the 
eigenf unctions e k satisfy 

e k G C (D), \e k (x,y)\ < C, 

\d x e k (x,y)\, \d y e k (x,y)\ < Cy/\X^\, 

for (x, y) G D, positive integer k, and some constant C > 0. For the square 
domain D = (0, 1) x (0, 1), these conditions are all satisfied. Then, according 
to Theorem 5.2.9 in the stochastic convolution Z(t) has a continuous 
version with values in L 2 (D). (Actually, in this case, Z(t) is in Cq(D), the 
Banach space of continuous functions satisfying the zero Dirichlet bound- 
ary condition on D). For this Wiener process W(t) in (||), the stochastic 
convolution Z(t) is 

oo . t 

Z{t) = Y^a k e k / e-W-Up^s) , t>0. (10) 
k=i Jo 

As shown in 0, for every initial condition o>(0) G L 2 (D), there exists a 
unique global mild solution uj(x,y,t) of the quasigeostrophic flow model (|3|). 
This solution is in C([0, T];L 2 (D)) for every T > 0. 

3 Existence of an Invariant Measure 

Now we consider invariant measure for the quasigeostrophic flow model (|3|). 
For the rest of the paper, we denote u(t; x) as the solution of the quasi- 
geostrophic flow model with initial condition (not the spatial point) x G H. 
We introduce the usual notations. The Markovian transition semigroup 

is 

(P t g)(x)=E[g(u(t; x)% 

for g G Bb(H), the space of bounded Borel measurable functions. Hereafter 
E is the expectation. The transition probability is 

p t (x,r)=P(w(t; x)er), 
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for x G H and T G 13(H), the cr-algebra of Borel sets in H. 
A probability measure \x on (H, 13(H)) is called invariant if 



gdjx = J Ptgdfj, 
for any t > and g G Bb(H), or, equivalently, 

P t (x,T)dfj J = n(F), 

H 

for any t > 0, x G # and T G -B(-ff). 

The existence of an invariant measure for the quasigeostrophic flow model 
(H) follows from a tightness or , equivalently, a compactness argument ET 



If the mean-square norm of the solution is bounded for all time t > and for 
all initial data, then by the Chebyshev inequality, the solution is bounded in 
probability, which further implies that the family of measures on (H, 13(H)) 



l -j P t (x,-)dt, T>1, 



T 

is tight for some x G H; see ||, page 89-90. Thus by Corollary 3.1.2 in J3|, 
there exists an invariant measure for the quasigeostrophic flow model (|3|). So 
in the rest of this section, we estimate the mean-square norm E||cj(t)|| 2 . 
We assume that 

POO 

/ \\S(r)VQ\\ 2 HS dr < +oo, (11) 
Jo 

where || ■ \\jjs is the Hilbert-Schmidt norm. We rewrite (|5|) as 

u(t) = Y(t) + Z(t) , (12) 

where ^ 

Y(t) = S(t)x+ [ S(t- s)F(uj(s))ds, 



o 



with initial data u>(0) = x, and Z(t) is the Ornstein-Uhlenbeck process in 
®- 

By Corollary 4.14 in 0, for any x G H, 

supE||Z(t)|| 2 = supE f ||5 , (r)v / Q||H S dr<+oo. (13) 

t>0 t>0 Jo 
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By or follow a Yosida approximation combined with L 2 -norm estimate 
as in Proposition 6.1.6 in [ffl, we have, for any x G H, 



Note that 



\ui 



supE||y 

t>0 



-<Y + Z,Y + Z > 



< +00. 



(14) 



< 
< 



||F|| 2 + 2 < Y,Z > 

||r|| 2 + 2||r|| ||z|| -\ 
2fi|y|| 2 + iizii 2 ). 



\z\ 



Thus, by (0) and flTJ, 



sup E,\\lo 

t>0 



< +00. 



(15) 



By the argument in the beginning of this section, there exists at least one 
invariant measure for the the quasigeostrophic flow model (0). We have the 
main result in this section. 

Theorem 1 Assume that / + °° ^S^^/Q^gdr < +00. Then there exists at 
least one invariant probability measure for the quasigeostrophic flow model 
(fjj in the space L 2 (D) of square-integrable vorticities. 



4 Uniqueness of an Invariant Measure 

Now we consider the uniqueness of invariant measure for the quasigeostrophic 
flow model (|3|). As we know in Chapter 4 in 01] , the uniqueness of invari- 
ant measure is a consequence of regularity of the transition semigroup P t , 
by the Doob's Theorem. Due to Khasminskii's Theorem, strong Feller and 
irreducibility properties imply the regularity. So we now try to prove the 
strong Feller and irreducibility properties for the transition semigroup Pt- 

Strong Feller property means that for every g(x) in Bf,(H), the space of 
bounded Borel measurable functions on H, Ptg{x) is in Cb(H), the space of 
bounded continuous functions on H. 

Irreducibility property means that for every Borel set in H, i.e., for every 
T in B(H), Pt(x, T) is positive for any x G H and t > 0. 



Strong Feller Property 
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We first consider strong Feller property. Note that (||, p. 119) 

Trace f S{r)QS*{r)dr = [ \\S(r)^Q\\ 2 HS dr. 
Jo Jo 

So the condition for the existence of invariant measures in Theorem |l|, i.e., 
Jo + °° ll^( r ) VQ\\ 2 Hs^ r < +°°; implies that the linear integral operator Q t : 

Q t x := I S(r)QS*(r)xdr, x G H, (16) 
Jo 

is of trace class for any t > 0. 
We further assume that 

i 

Image S(t) C Image Ql . (17) 

Then follow a similar argument as in the proofs of Theorem 7.2.4 in |IJ and of 
Theorem 3.1 in |)TU|, we conclude that P t , t > 0, is a strong Feller semigroup. 



Irreducibility Property 

Now we consider irreducibility property. We further assume that the 
covariance operator Q is one-to-one (or injective), i.e., the kernel ker Q = {0}. 
Then, as in the proof of Theorem 7.4.2 in || and of Theorem 3.1 in |10| , P t , 
t > 0, is irreducible. 



Thus, with the strong Feller and irreducibility properties proved above, 
using Doob's Theorem 4.2.1 in Hj, there exists a unique invariant measure 
on (H,B(H)), and all other transition probability measures Pt(x, ■), x £ H, 
approach this unique invariant measure // as time goes to infinity. 

Therefore, we have the following main theorem in this section. 

Theorem 2 Assume that 

(i) f +Oa \\S{r)y/Q\\ 2 HS dr < +oo ; 

i 

(ii) Image S(t) C Image Ql , where Q t is defined in fildj), and 
(Hi) The covariance operator Q : L 2 (D) — > L 2 (D) is one-to-one. 
Then 
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(A) There exists a unique invariant probability measure \i for the quasi- 
geostrophic flow system (|3[) in the space L 2 (D) of square-integrable vorticities; 

(B) Moreover, for any uj £ L 2 (D), the transition probability measures 
Pt{w, ■) approach the unique invariant probability measure \i. Namely, for 
any T £ B{H), 

lim p t ( Wj r) = ^(r) ; 

t— »+oo 

and 

(C) Quasigeostrophic flow system ^) is ergodic, namely, 
lim / g(u(t))dt = / gdjj,, P — a.s. 



L- 



for all solution u(t) with initial date in L 2 (D) and all Borel measurable func- 
tion g : L 2 (D) — > R such that J L i, D ^ ||p||<^/-i < 00 ■ 

The ergodicity in Part (C) above is a consequence of the uniqueness of 
the invariant measure jx] see Theorem 3.2.6 in H. 



5 Summary 

In this paper, we have studied ergodicity of large scale quasigeostrophic flows 
under random wind forcing. We have shown that the quasigeostrophic flows 
are ergodic under suitable conditions on the random forcing and on the fluid 
domain, and under no restrictions on viscosity, Ekman constant or Coriolis 
parameter. When these conditions are satisfied, then for any observable 
of the quasigeostrophic flows, its time average approximates the statistical 
ensemble average, as long as the time interval is sufficiently long. 

There is recent work on random dynamical attractors for the quasi- 
geostrophic flow model by Duan et al. |J. A consequence of that work 
implies that, when viscosity is sufficiently large and when the trace of the 
covariance operator for the Wiener process is sufficiently small, then all quasi- 
geostrophic motions approach a point random attractor exponentially fast as 
time goes to infinity. This is a very rare case. This point random attractor 
corresponds to a unique invariant Dirac measure, i.e., the supporting point 
of the Dirac measure is the global (point) attractor, and thus under these 
conditions, quasigeostrophic flows are also ergodic. These conditions are dif- 
ferent from the ergodic conditions in the current paper. For example, in the 
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current paper, we do not impose any condition on viscosity, or on the size of 
the trace of the covariance operator for the Wiener process. 
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